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Abstract. We describe a new Maple package for treating boundary 
problems for linear ordinary differential equations, allowing two-/multi- 
point as well as Stieltjes boundary conditions. For expressing differential 
operators, boundary conditions, and Green's operators, we employ the 
algebra of integro-differential operators. The operations implemented for 
regular boundary problems include computing Green's operators as well 
as composing and factoring boundary problems. Our symbolic approach 
to singular boundary problems is new; it provides algorithms for com- 
puting compatibility conditions and generalized Green's operators. 
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1 Introduction 

Although boundary problems clearly play an important role in applications and 
in Scientific Computing, there is no systematic support for solving them symbol- 
ically in current computer algebra systems. In this paper, we describe a Maple 
package with algorithms for regular as well as singular boundary problems for lin- 
ear ordinary differential equations (LODEs). While a first version of the package 
with functions for regular boundary problems was presented in [T], the methods 
and the implementation for singular problems are new. A prototype implemen- 
tation for regular boundary problems in the TH30REMV system was described 
in [2] as part of a general symbolic framework for boundary problems, including 
also some first steps towards linear partial differential equations (LPDEs). 
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In Section [5J we recall the algebra of integro-differential operators provid- 
ing the algebraic structure for computing with boundary problems. We describe 
its implementation in Maple, where we use a normal form approach in con- 
trast to In Section [3l we outline our symbolic approach for solving boundary 
problems. For an analytic treatment of boundary problems for LODEs, see for 
example f3l4] or [5] for further applications. The functions we present include 
the computation of Green's operators and Green's functions as well as the fac- 
torization of boundary problems. 

We introduce generalized boundary problems in Section 2] and develop an 
algorithm for computing generalized Green's operators. The main step of the 
algorithm is to determine compatibility conditions for arbitrary boundary prob- 
lems in an algebraic setting; the special case of two-point boundary problems of 
second order is discussed in [6, Lecture 34]. For singular boundary problems and 
generalized or modified Green's functions in Analysis, we refer for example to [1] 
and [7], and in the context of generahzed inverses to [51 Sect. 9.4], [5], and [TOl 
Sect. H]. 

The Maple package IntDiffOp is available with an example worksheet at 
http : //www . rise . jku . at /people/ akorpora/ index . html 



2 Integro-Differential Operators 

We first recall the definition of integro-differential algebras and operators, see 
PT| and [12 for further details. For the similar notion of differential Rota-Baxter 
algebras, we refer to [13]. As a motivating example, consider the algebra F — 
C°°(]R) with the usual derivation and the integral operator /:/>-> /(^) for 
a fixed a G K. The essential algebraic identities satisfied by the derivation and the 
integral operator are the Leibniz rule, the Fundamental Theorem of Calculus, 
and Integration by Parts. Note also that /(a) = / — J/', so the evaluation 
Ea : / I— ^ fio) at the initialization point a of the integral can also be expressed 
in terms of the derivation and integral. 

We call (7^, (9, J) an integro-differential algebra if {J',d) is a commutative 
differential algebra over a commutative ring K and J is a if-lincar right inverse 
(section) oi d = ' , meaning (J/)' — f, such that the differential Baxter axiom 

(J/')(J.9') + J(/.9)'-(J/')3 + /(j3') 

holds. We call e = 1 — Jo9 the evaluation of J'. We say that an integro-differential 
algebra over a field K is ordinary if Ker(9) = K. For an ordinary integro- 
differential algebra, the evaluation can be interpreted as a multiplicative linear 
functional (character) e: ^ K. This allows treating initial value problems, 
but for doing boundary problems we need additional characters ip: ^ K {in 
the above example, evaluations Ec: / i— > /(c) at various points c £ K). 

Let {T, d, J) be an ordinary integro-differential algebra over a field K and 
let C J^* be a set of multiplicative linear functionals ip: J- K including e. 
The integro-differential operators J-<p[d, are defined in [IT] as the if-algebra 
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Table 1. Rewrite Rules for Integro-Differential Operators 



generated by the symbols d and J, the "functions" f € J- and the "functionals" 
(f G "P, modulo the Noetherian and confluent rewrite system of Tabic [1] 

The representation of integro-differential operators in our Maple implemen- 
tation is based on the fact that every integro-differential operator has a unique 
normal form as a sum of a differential, integral, and boundary operator. The 
normal forms of differential operators are as usual /i^N integral operators 
can be written uniquely (up to bilinearity) as sums of terms of the form f Jg, 
and the normal forms of boundary operators are given by 

/^v^^' + 5j.vV'J/»i,¥>) , (1) 

with only finitely nonzero summands. Stieltjes boundary conditions are boundary 
operators where fi^^p = a^p^i G K and gj^tp = 1. They act on !F as linear functionals 
in the dual space . See [2] for Stieltjes boundary conditions in Analysis. 

From Table [1] formulas can be derived for expressing the product of integro- 
differential operators directly in terms of normal forms; see |15| for the case 
# = {e}. Implementing these formulas leads to faster computations since we 
need not reduce in each step. In our package, we use for the underlying "integro- 
differential algebra" all the smooth functions in one variable representable in 
Maple, together with the usual derivation and the integral operator J = Jg, 
both computed by Maple . We take as characters # = {Ec | c G R}. 

We created data types for the different kinds of operator, representing integro- 
differential operators as triples iNTDiFFDp(a, 6, c), where a is a differential oper- 
ator, b an integral operator and c a boundary operator. Differential operators are 
represented as hsts diffop(/o, /i, . . .) and integral operators as lists of pairs of 
the form intop(intterm(/i, gi), intterm(/2, 32), • • •)• In order to have a unique 
representation for integral operators, one would need a basis of the underlying 
integro-differential algebra and use only basis elements for the gi. In our imple- 
mentation, we use the following heuristic approach: We split sums in the gi and 
move scalar factors to the coefficients fi. 

Due to (dl), a boundary operator bdunddp contains a list of evaluations at dif- 
ferent points. Each evaluation evop is a triple containing the evaluation point, the 
local part ^ fi,ipfd^ and the global part ^ gj^^ip^hj^^. Hence we use the expres- 
sion boundop(evdp(c, evdiffdp(/o, . . .), evintop(evintterm((7i, hi), ...),...) for 
the representation of boundary operators. 

In the following example, we first enter some operators of different types. For 
displaying the operators, we use D for d, A for J and E[c] for the evaluation e^. 
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> T 


:= DIFFDP(0,0,1) ; 




r := 


> G 


:= INT0P(INTTERM(1,1)) ; 




G := A 


> B 


:= B0UND0P(EV0P(1, EVDIFFDP(l), EVINT0P(EVINTTERM(1 , 1) )) ) ; 




B := E[l] + ((E[ll) . A) 



Now we show how to add and multiply integro-differential operators and how 
to apply them to a function f G T. 



> 


ApplyOperator(G, f(x)); 






> 


MultiplyOperator(G,G) ; 




(x . A) - {A . x) 


> 


MultiplyOperator(T,G,G) ; 

1 


> 


S := AddOperatord, G, B) ; 




S ~ -D^ + A + E[l] + ((E[l]) . A) 


> 


ApplyOperator(S, f(x)); 




|^/(^) + Jo/(^)rf^ + /(I) + Ilfi^)dx 
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In this section, we demonstrate how to compute with regular boundary problems 
in our MAPLE package. For an integro-differential algebra a boundary problem 
is given by a monic differential operator T = 9" + c„_i9"^^ + ■ • ■ + cid + cq 
and boundary conditions . . . ,/?m- Given a forcing function / G J^, we want 
to find u ^ J- such that 



Tu = /, 




PiU = ■ 


■ = PnU = 0. 



(2) 



A boundary problem is called regular if for each f (z J- there is exactly one u ^ T 
satisfying We want to solve a boundary problem not only for a fixed / but 
to compute the Green's operator mapping each forcing function / to its unique 
solution u. In other words, we solve a whole family of inhomogeneous differential 
equations, parameterized by a "symbolic" right-hand side /. We restrict our- 
selves to homogeneous conditions because the general solution is then obtained 
by adding a particular solution satisfying the inhomogeneous conditions. 

For convenience, we shortly recall the abstract linear algebra setting for 
boundary problems over a vector space T as described in [K]. For U < T 
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we define the orthogonal as IJ-^ — {13 e T* : /3(ii) = for all u e [/} < J"*. 
Similarly, for ^ < J"*, we define = {v e T : l3{v) = for all ^ e B} < T. A 
subspace U (resp. B) is orthogonally closed ii U = If-'--'- (resp. B = B^^). Every 
subspace U < J- is orthogonally closed and every finite dimensional subspace 
B < J-'* is orthogonally closed. For a linear map T: T ^ Q between vector 
spaces, the transpose map T* : Q* ^ T* is defined by 7 i— >■ 7 o T. The image of 
an orthogonally closed space under the transpose map is orthogonally closed. 

A boundary problem is given by a pair (T,B), where T is a surjective linear 
map and B < J^* is an orthogonally closed subspace of the dual space. We call 
It e J-" a solution of (T, B) for a given / G if Tu = f and u e B^. A boundary 
problem is regular if for each / there exists a unique solution u. The Green's 
operator of a regular problem maps each / to its unique solution u. We also 
write {T,B)~^ for the Green's operator. A boundary problem is regular iff B'^ 
is a complement of Ker T so that T = Ker T + B^ as a direct sum. 

For F = C°°[a, 6], a monic differential operator T is always surjective and 
dim Ker T — n < 00. Moreover, variation of constants can be used to compute 
a distinguished right inverse: If T has order n and ui, . . . , u„ is a fundamental 
system for it, the fundamental right inverse is given by 



where d is the determinant of the Wronskian matrix W for (ui, . . . , it„) and di 
the determinant of the matrix Wi obtained from W by replacing the i-th column 
by the n-th unit vector. Equation ^ is valid in arbitrary integro-difFerential al- 
gebras provided the n-th order operator T has a fundamental system . . . , it„) 
with invertible Wronskian matrix; see or |12j . This will be assumed from now 
on, together with the condition dimB < 00 appropriate for LODEs. 

Regularity of a boundary problem (T, B) can be tested algorithmically as 
follows. If (ui, . . . , Un) is a basis for KerT and . . . , for B, we have a 
regular problem iff the evaluation matrix 



is regular; see [ini Cor. A. 17] or [T71 p. 184] for the special case of two-point 
boundary conditions. Of course this implies m = n, but we will consider more 
general types of boundary problems in Section 2] where this is no longer the case. 
It will also be convenient to use the notation ^ for arbitrary iti, . . . ,m„ G T 
and boundary conditions . . . , 

The algorithm for computing the Green's operator is described in detail in 
[11] : see also [2]. The main steps consist in computing the fundamental right 
inverse £ T[d, J] from a given fundamental system as in ^ and the projector 
P G T[d, J] onto KerT along B-^. Then the Green's operator is then computed 




(3) 





as G = (1 - P)T*. 



6 A. Korporal, G. Regensburger, M. Rosenkranz 



For a boundary problem we need to enter a monic differential operator T 
and a list of boundary conditions (51, . . . , hm) as described in Section [5] in the 
form Bp(r, Bc(61, . . . , bm)). We use the Maple function dsolve for computing a 
fundamental system of T . As an example, we compute the Green's operator for 
the simplest two-point boundary problem u" ~ /, u(0) = u{l) = 0. From the 
Green's operator for two-point boundary problems, we can extract the Green's 
function [18 , which is usually used in Analysis to represent the Green's operator. 

> T := DIFFQP(0,0,1) : 

> bl := BOUNDOP(EVOP(0, EVDIFFOP(l), EVINTQPO)): 

> b2 := BQUND0P(EV0P(1, EVDIFFDP(l), EVINTQPO)): 

> Bp := BP(T, BC(bl, b2)) ; 

Bp ■- BP(D^ BC(E[0], E[l])) 

> IsRegular(Bp) ; 

true 

> GreensOperator (Bp) ; 

(a; . A) - (A . a;) - {{x E[l]) . A) -I- {{x E[l]) . A . a;) 

> GreensFunctionC/i) ; 

+ x^ ^ and ^ <= x and x <= 1 
—X + xS, X and x <— ^ and ^ <= 1 

For simplifying boundary problems, we can apply factorizations into lower 
order problems along given factorizations of the differential operators. Further 
details and proofs of the following results can be found in [TB] and [TT]. The 
composition of two boundary problems (Ti, Bi) and (T2, B2) is defined as 

(ri,Bi) o {T2,B2) - {TiT2,t;{Bi)+B2). (5) 

The composition (Ti, Bi) o (T2, /B2) of two regular boundary problems is regular 
with Green's operator 

{{Ti,Bi)o{T2,B2))-^ = {T2,B2)-HTi,Bi)-\ (6) 

Given a regular boundary problem {T,B)^ every factorization T — T1T2 can be 
lifted to a factorization {T,B) = {Ti,Bi)o (T2,B2), where (Ti,Bi) and (T2,S2) 
are regular and B2 < B. For factorizing a differential operator, we use the func- 
tion DFactor in the Maple package DEtools. As an easy example, we show 
how to factor the boundary problem from above; more examples for solving and 
factoring boundary problems can be found in our example worksheet. 



> Bp := BP(T, BC(bl, b2)) ; 




Bp ■- BP(D^ BC(E[0], 


E[l])) 


> fl, f2 := FactorBoundaryProblem(Bp) ; 




/I, /2 := BP(D, BC(E[1] . A)), 


BP(D, BC(E[0])) 
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4 Singular Boundary Problems 

For illustrating the main issues with singular boundary problems, we consider 
the boundary problem 



u" ^ /, 

m'(0) = u'{l) = 0; 



(7) 



see for example [U Page 215] or [THl Section 3.5] from a Symbolic Computation 
perspective. This problem is singular since it is not solvable for all / G J^. It can 
easily be seen that if u" = /, then / has to fulfill the compatibility condition 
u'{l) = Jq /(^) dS, = 0. Moreover, uniqueness fails as well: If a solution u £ T 
exists, then also u + c solves the problem for all c G M. 

Our goal here is to generalize the symbolic approach of the previous section 
to problems of the kind Since we want to compute generalized Green's op- 
erators, we cannot give up uniqueness of solutions — but we no longer require 
existence. Of course, uniqueness of solutions can always be achieved by impos- 
ing additional boundary conditions. On the other hand, adding too many condi- 
tions introduces new compatibility conditions, which we want to avoid (see after 
Lemma [1] for the precise statement). For the boundary problem ([7]), we can add 
for example the condition u(l) = and consider the problem 



u'{0) = u'{l) = u(l) = 0. 



(8) 



This does not introduce any new compatibility conditions as we will see later 

(see before Lemma [5]). 

A boundary problem has at most one solution for each forcing function / iff 
n KerT = {0}. We see that for © we have 6^ n KerT = R while in © the 

intersection is {0}. The regularity test for boundary problems in terms of the 

evaluation matrix Q can be generalized from the setting in Section [3J 

Lemma 1. Let U — [ui, . . . , m„] < T and B = . . . , /3„i] < T* with f3i and 
Uj linearly independent. Then U D = {0} ijf the evaluation matrix f3{u) has 
full column rank. 

Proof. Let bj denote the columns of /3(u). The evaluation matrix has deficient 
column rank iff there exists a linear combination X]j=i '^j^j ~ ^ with at least one 
Xj 0. This is the case iff there exist a nonzero u — X^jLi ^j^j ^ U O B^. □ 

As mentioned for the example ([5]), singular boundary problems typically 
impose compatibility conditions on the admissible forcing functions. We can now 
make this precise: Clearly, a function / is admissible iff it is of the form Tu 
for a function u that satisfies the boundary conditions from B, so the space of 
admissible functions is T{B^). The compatibility conditions provide an implicit 
description of this space, comprising all those linear functionals that annihilate 
T{B'^). In other words, the compatibility conditions are the subspace T{B'^)^ 
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of T* . This also makes precise what we mean by adding boundary conditions 
without imposing additional compatibility conditions: We enlarge to B so as 
to ensure n KerT = {0} despite retaining T{B^) = T{B^). 

For tackling the problem of existence, we modify the forcing function. In 
the example this looks as follows: Since a solution exists only for forcing 
functions that fulfill f{Cj d.$, — 0, we consider the problem 



u'{0) = = u{\) = 0, 



(9) 



which now always has a unique solution. For those / that fulfill the compatibility 
condition, problem (|5]) remains unchanged. 

The general idea is that we project an arbitrary forcing function into the space 
of admissible functions. But this involves choosing those "exceptional functions" 
that we want to filter out. Even in the simple example ([5]), we might as well 

project / to / — 52^X0/(0 '^C instead of / — Jq/(0 In the second case, we have 
filtered out the constant functions, in the first case the linear-homogeneous ones. 
The space £ of exceptional functions can be any complement of the space T(i3^) 
of admissible functions, like £ = [1] or £ = \x\ in this example. 

Definition 1. A generalized boundary problem is given by a triple {T,B,£), 
where (T, B) is a boundary problem and £ <T . A generalized boundary problem 
is called regular if 

S^nKerT={0} and F = T{B^)+£. 

The generalized Green's operator maps each forcing function f to the unique 
solution of the boundary problem 



Tu = Q/, 
P^u =...-- 



PrnU = 0, 



where B = . . . , /?,„] and Q is the projector onto T[B^) along £. We also 
write {T, B, £)^^ for the Green's operator. 

If {T,B,£) is regular, the restriction T|gj_ : B-^ T{B'^) is bijective. So the 
generalized Green's operator is given by 

1-1 ^ 



G = T\ 



(10) 



We begin with computing the projector Q. For this we derive first an explicit 
description of the space of compatibility conditions. 

Proposition 1. Let {T,B,£) be a generalized boundary problem and let G be 
any right inverse of T . Then we have 

T{B^)^ = G*{Bn{KcrT)^). (11) 

Moreover, dimr(i3^)^ = dim£ for any complement £ with J-' — T{B^) + £. 
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Proof. With J£j Prop. A.6], we see that T{B^)^ = (T*)-i(i3). Since T is sur- 
jective, T* is injective, and for any right inverse G of T, G* is a left inverse of 
T*. Hence {T*)-\B) = G*{BnlmT*) by [H Prop. A.13]. Again by [H Prop. 
A.6], we have ImT* = (KerT)-^, and hence T{B-^)-^ = G* (B H (Ker T)^) . 

Since dimS < oo, by the first statement also dimr(;B-'-)-'- < oo. But T{B^) 
is orthogonally closed; see for example Section A.l]. Therefore we obtain 

dimr(B^)^ = codimr(i3-L)^^ = codimT(B^), 

and the statement follows immediately from [16, Prop. A. 14]. □ 

Note that s — dimf = codimr(;B^) counts the number of (linearly inde- 
pendent) compatibility conditions. Equation ([TT|) is the key for an algorithmic 
description of the projector Q onto T{B^) along £. The space E is given as part 
of the problem description, and it can be specified by a basis (wi, . . . , Ws). Since 
the other space T{B^) has finite codimension s, it can be specified in terms of s 
linearly independent compatibility conditions, and Equation (|lll) can be used 
to compute these in terms of T and B. For that we just have to determine a 
basis of B n (KerT)^ and then apply any right inverse G of T, for example the 
fundamental right inverse defined in Section [3l 

For determining a basis of ;B n (KerT)-*- we first compute the kernel of the 
transpose of the evaluation matrix /?(u), where (wi, . . . ,it„) is any basis of KerT 
and (/?!,..., Pro) any basis of B.li w — (wi, . . . , WmY & Ker (3{uy, then 

m 

«;*(/3i,...,/3™)* eSn(Kerr)^, 
t=i 

hence a basis of B n (KerT)-'- can be obtained by computing the products 
(u*i(/3i, . . . ,;0m)*, . . ■ ,Wfc(/3i, ■ . • ,/3m)*), where (ui, . ..,«&) is a basis of Ker /?(«)*. 

Using Proposition [U we can now verify that the compatibility conditions of 
the boundary problems ^ and ([8]) are the same. In both cases we have T = d^, 
so we can choose the fundamental right inverse J J = xj — Ja; and (l,x) as a 
basis of Ker T. The evaluation matrices are given by 

/3(w) = (^0 J ) and (3{u) 

In the first case, a basis of /3(u)* is given by ((—1, 1)*), hence (Ei9 — EqO) is a 
basis of ;6n(KerT)-'-. In the second case, a basis of /3(u)* is given by ((—1, 1, 0)*) 
and the basis of S n (KcrT)^ is again (Ei9 — Eo9). Multiplying this basis by 
the right inverse of T, we get as a basis for the compatibility conditions 

{Eld - Eod) ■ {xj - Ja;) = Ei{xd + 1)J - Eo{xd + 1)J - Eid^x + Eodjx 
= Eix + El J - Eox - EoJ - Eix + Eax = Ei J = J^, 

which agrees with our heuristic considerations after ([7]). 
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We can now compute the projector Q just as the kernel projector P for 
standard boundary problems (mentioned in Section If . . . , Kg) is a basis 
for the compatibility conditions T{B'^)'^ and {wi, . . . ,Ws) a basis for £, then 
the corresponding evaluation matrix k{'w) is regular by Lemma [U which can be 
applied to J" = T{B^) +£ = T{B^)^^ + £ since T{B^) is orthogonally closed. 
Hence we can compute the projector Q onto T{B^) along £ as 

s 

(9 = 1 - y^^wjkj, 

where . . . , hsY — n{w)~^ ■ (ki, . . . , Ks)*; see for example [16l Lemma A.l]. 

The final step for computing the generalized Green's operator pUj) is to find 
the inverse function T\j^]_. In the regular case, we started with an arbitrary right 
inverse of T and multiplied with a projection onto B^ along KerT. But this step 
cannot be generalized to our setting. Our approach is to embed the generalized 
problem into a standard one in the following sense. 

First note that the evaluation matrix of a regular generalized boundary prob- 
lem has full column rank by Lemma [TJ so it has a left inverse. 

Lemma 2. Let (T, B, £) be a regular generalized boundary problem. Let P{u)~ be 
a left inverse of j3{u) and (/3i, . . . , /?„)* — /3{u)~ (Pi, . . . , PmY ■ Then the boundary 
problem (T, B) is regular, where B < B is spanned by /3i , . . . , /3„ . 

The proof of the statement is obvious, since the evaluation matrix /3(u) is 
given by j3{u)~ (3{u) = 1„. Hence the problem (T, is regular. In our package, 
we always choose the Moore- Penrose pseudoinverse as a left inverse (5{u)~ of 
the evaluation matrix P[u). The generalized boundary problem for example 
embeds into the standard boundary problem 



u" = /, 

it'(O) V u'(l) - 2 u(l) = u'(0) u'(l) = 0. 

The Green's operator for this regular problem according to Section [3] is given 
by a; J — Ja; — ^(x -I- 1) -I- J^a;. The next proposition tells us how to compute the 
generalized Green' s operator from it. 

Proposition 2. Let {T,B,£) be a regular generalized boundary problem and let 
{T,B) be a regular boundary problem with B < B. Then 

{T,B,£)-' = {T,B)-'Q, 

where Q is the projector onto T{B^) along £. 

Proof. Since < ;B, we have B^ < B-^. Hence the maps T\^i and G = (T, B^^ 
coincide on Since Tjgi : B^ — s> T{B^) is a bijcction, we can compute the 
restriction T\~^\ by first applying a projector onto T{B^) and then G. Hence 
T\-i = GQ, where Q is again the projection onto T{B^) along £. Hence the 
generalized Green's operator is given hy G — T\„^Q = GQ^ = GQ. □ 



(12) 
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Applying the previous proposition to Example ([S]) leads to the generahzed 
Green's operator a;J — Jx — ^(a;^ + l)Jo + Jo^- ^ more involved example 
illustrating the Maple functions in our package, we refer to the Appendix. 

5 Outlook 

We are currently investigating in how far the composition of boundary prob- 
lems ([5]) can be extended to generalized boundary problems such that an analog 
of the "reverse order law" ^ holds. We can see in the example below that 
for such a generalization, we also have to modify the second component with 
the boundary conditions. The question under which conditions a reverse order 
law holds for different classes of generalized inverses — not necessarily related 
to integro-differential operators — is extensively studied in the literature, see for 
example [12] and the references therein. 

The search for generalized composition laws is intimately connected with the 
question of "embedding" a singular boundary problem into a regular problem 
of higher order. For example in [T8] j the Green's operator G of the generalized 
boundary problem (9^, [Eo9, Ei9, Jp], [1]) can be factored as G = G od where G 
is the standard Green's operator of the boundary problem (9^, [Eo9, Ei9, j^]). 
Hence G = G o and, assuming © for the composition, also 

{d^ [Eod, Ei9, Jj], [0]) = {d, [Eo], [0]) o id', [Eod, E,d, jj], [1]), (13) 

since Jq is the Green's operator of the boundary problem {d, [Eq]). The singular 
second-order problem is thus embedded into a regular third-order one. 

Multi-point boundary problems can also be treated by our method, yielding a 
suitable Green's operator just as in the classical two-point setting. Generalizing 
the extraction procedure for Green's functions is future work, see [50] for an 
analytic description of Green's functions for multi-point boundary problems. 

Going from LODEs to LPDEs, more drastic changes are necessary since ge- 
ometry enters the picture. For example, the Green's operator of the inhomoge- 
neous wave equation u^x—utt = f with homogeneous Dirichlet data on the x-axis 
integrates / over a certain triangle whose tip is at {x, t). In terms of the operator 
algebra, this means one must incorporate the chain and substitution rule along 
with explicit operators encoding change of variables. A first approach along these 
lines, for the very simple case of linear coordinate changes, was presented in [2] 
and is currently being refined. Studying singular boundary problems for LPDEs 
from a symbolic point of view is also very interesting; see for example |21) for a 
Grobner bases approach to compute the (hierarchy of) compatibility conditions 
for elliptic boundary problems. It would be tempting to combine the tools of 
involutive systems used there with the setting of operator rings used here. 
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A Example 

Now wc will give a detailed example for computations with generalized boundary 
problems. We introduced a new datatype gbp(T, bc(61, . . . , bm), es(/1, . . . , fk)), 
where T and (61, ... , bm) again are a differential operator and boundary con- 
ditions and (/I, . . . , /to) is a basis of the exceptional space. Wc added the new 
procedures Compatibility Conditions, IsComplement and Projector, which will 
be explained later and extended the procedures Greens Operator and IsRegular. 
The first one now also computes the Green's Operator for a generalized bound- 
ary problem and the second one tests the condition KerT fl B-^ = {0} also for 
generalized boundary problems. 

We consider the more complicated example 



u"" + u" = f 

u'(0) = u"(0) = u"{Tr) = u'"{0) = u'"{ir) = 0. 



(14) 



We enter the boundary problem stated above and compute a fundamental system 
for the differential operator T = + . 



T := DIFFOP(0, 0, 1, 0, 1) : 
b[l] := BOUNDOP(EVDP(0, EVDIFFOP(0, 1), 
b[2] := BOUNDOP(EVOP(0, EVDIFFOP(0, 0, 
b[3] := BOUNDOP(EVOP(0, EVDIFFOP(0, 0, 
b[4] := BOUNDOP(EVDP(Pi, EVDIFFOP(0, 0, 
b[5] := BOUNDOP(EVDP(Pi, EVDIFFOP(0, 0, 
Bp := BP(T, BC(b[l],b[2],b[3],b[4],b[5])): 
fs := FuiidcimentalSystem(T) ; 



EVINTOPO)) : 
1) , EVINTOPO)) : 
0, 1), EVINTOPO)): 
1), EVINTOPO)): 
0, 1), EVINTOPO)): 



[x, sin(a;), cos(a;), 1] 



Now we add another boundary condition b[6] in order to achieve uniqueness of 
solutions. This can be checked by considering the column rank of the evaluation 
matrix. We further verify that the compatibility conditions of both problems are 

the same. 



> b[6] := BOUNDOP(EVOP(Pi, ZEROEDOP, EVINTOP (EVINTTERM ( 1 , 1 ) ) ) ) : 

> BpA := BP(T, BC(b[l],b[2],b[3].b[4],b[5],b[6])): 

> IsRegular (BpA) ; 

true 

> CompatibilityConditionsCBp) ; 

BC((E[Pi] . A . (sin(a;)), (E[Pi]) . A . (cos(a;))) 

> CompatibilityConditions(BpA) ; 

BC((E[Pi] . A . (sin(a;)), (E[Pi]) . A . (cos(a;))) 
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Now we enter a generalized boundary problem and check that our choice [1, x] 
as exceptional space is a complement of T{B-^). Then we compute the projector 
Q onto T(B^) and the Green's operator for the generalized boundary problem 

(T, [b[l],b[2],b[3],b[4],b[5],b[6]],[l,x]), 

> gBp := GBP(T, BC(b[l],b[2],b[3],b[4],b[5],b[6]), ES(l,x)): 

> IsComplement (gBp) ; 

true 

> Q := Projector (gBp) : 

Q:=l- l((E[Pi]) . A . (sin(x))) + 1^ (-^ + I) . (E[Pi]) . A . (cos(x) 

> G := GreensOperator(gBp) : 



Finally we verify that the Green's operator G fulfills the equation TG = Q and 
the six boundary conditions. 

> simplify (SubtractOperator(MultiplyOperator(T, G) , Q)) 



> seq(simplify(MultiplyOperator(b[i] , G)), i=1..6); 

0, 0, 0, 0, 0, 



